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Abstract. We study the bidimensional voter model on a square lattice with 
numerical simulations. We demonstrate that the evolution takes place in two 
distinct dynamic regimes; a first approach towards critical site percolation and 
a further approach towards full consensus. We calculate the time-dependence 
of the two growing lengths finding that they are both algebraic though with 
different exponents (apart from possible logarithmic corrections). We analyse the 
morphology and statistics of clusters of voters with the same opinion. We compare 
these results to the ones for curvature driven two-dimensional coarsening. 
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1. Introduction 

Purely dynamical stochastic models are used to describe problems beyond physics 
such as opinion formation pQ and population genetics [2], and treat issues in ecology, 
linguistics, etc. In the former context, questions on the spatial spreading of opinions 
are posed in terms of coarsening or segregation, just as in physical materials. 

The voter model laiaE] is one such purely dynamical stochastic model, used 
to describe the kinetics of catalytic reactions and as a prototype model of 

opinion dynamics P HU]. In its simplest realisation a bi-valued opinion variable. 
Si = ±1, is assigned to each site on a lattice or graph with some procedure that 
determines the initial conditions. Typically, the initial state is taken to be unbiased, 
with equal number of one and the other state. The dynamic rule is straightforward: 
at each time step, a variable chosen at random adopts the opinion of a randomly- 
chosen neighbour. These moves mimic the influence of the neighbourhood on the 
individual opinion. The probability of the chosen spin to flip in a time step is 
simply given by the fraction of neighbours with opposite orientation. The model 
is parameter free and invariant under global inversion of the spins, that is to say, 
Z 2 symmetric. As a site surrounded by others sharing the same opinion cannot 
fluctuate, there is no bulk noise and the dynamics are uniquely driven by interfacial 
noise. In some papers the model is dehned in terms of a site-occupation variable 
instead of a spin. 

Mathematicians, more precisely probabilists, solve this model by using a 
mapping to random walk theory [M- Physicists, instead, treat it within the master 
equation formalism. Once written in this form, one reckons that the transition 
probabilities do not satisfy detailed balance and, therefore, the model is essentially 
out of equilibrium. Even though there is no asymptotic thermal state, the dynamics 
can be solved as Glauber did for the stochastic Ising chain since the equations for 
the correlations of different order decouple PEI HD- 

The voter model’s evolution shows spatial clustering of similar opinions. It 
approaches one of the two absorbing states with complete consensus via a coarsening 
process in d < 2. It may also approach consensus in hnite size d > 2 systems but only 
because of a large random fluctuation with some small probability (that vanishes 
in the inhnite size limit). Otherwise, an inhnite family of disordered steady states 
exist in d > 2 mm The coarsening process in d = 2 is very different from the 
curvature driven one, as can be appreciated in Fig. where a series of snapshots 
of the spin conhguration at increasing times are shown, proving that the long-term 
dynamics are not determined by symmetry properties alone. It is also different from 
critical relaxation, specially because of the lack of bulk fluctuations, see also Fig. 
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The absence of bulk noise and surface tension entail important differences with 
respect to curvature driven phase ordering kinetics [niElDEI. In the voter model, 
regions of one opinion can only be penetrated by the other at the boundary. Besides, 
a large bubble consisting of voters of the same opinion does not shrink as in curvature 
driven processes. It slowly disintegrates as its boundary roughens diffusively to reach 
a typical width of the order of the initial radius [mini, while the radius of the 
droplet remains statistically constant (the radially-averaged magnetisation prohles 
have a stationary middle point). 

Coarsening processes usually conform to the dynamic scaling hypothesis nsna 
[T5] . This assumption states that if there is a single growing length in the process, say 
i{t), the statistical properties of the system are self-similar with respect to it. Under 
this assumption the space-time correlation is independent of time when distance is 
rescaled by i{t). In the voter model the evolution of a random initial condition shows 
the growth of ordered spatial regions. However, the exact asymptotic solution of 
an infinite system in d = 2 exhibits logarithmic violations of the standard scaling 
forms [13]. Although a characteristic length i(t) ~ y/i can be identihed, the density 
of interfaces decays as 1/ In f and the scaling function for the space-time correlation 
function C{r, t) involves an additional logarithmically decaying factor [TT] (somehow 
similarly to the critical dynamic scaling [T8l fT9] though with a logarithm instead of 
an algebraic correction). 

The goal of this work is to characterise the coarsening process in the 
bidimensional voter model with large but finite linear size by studying, in detail, the 
geometric and statistical properties of the dynamic pattern of domains. Following 
the analysis in [20] l2T] l22] 123] we will demonstrate that the system evolves in two time 
regimes: a pre-asymptotic approach to critical percolation and an ultimate approach 
to full consensus. With the aim of identifying and distinguishing the growing length 
in each of these regimes, we compute special time-dependent observables such as the 
number density of domains with a given area, or the number density of interfaces 
with a given length. As the characteristic length associated to the approach to 
percolation, that we call ip{t), grows quite slowly in time we are able to analyse this 
dynamic regime in detail (contrary to the what happens in the Ising model, where 
ip{t) is such a fast growing function of time that in practice critical percolation is 
reached too quickly to allow for a careful study of this dynamic regime). 

The paper is organised as follows. In Sec. we introduce the model and we 
summarise the time-dependence of several observables that were derived analytically 
in the past for inhnite size systems. In Sec. we present our numerical results. We 
discuss the violation of the scaling hypothesis as observed in the time-dependence of 
the density of interfaces, persistence probability, two-time correlation function and 
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space-time correlation function of an infinite size system. We then show our novel 
results on geometric properties of the largest cluster, number density of domain areas 
and interface lengths in hnite size systems. We relate their properties to the fractal 
properties of these objects. More importantly, this analysis allows us to demonstrate 
the existence of the two dynamic regimes evoked in the previous paragraph: a hrst 
approach to critical site percolation and the further evolution towards complete 
consensus in a longer time-scale. We end the paper with a concluding section. 


2. Analytical results 


The dehnition of the voter model is extremely simple. Each node i of a graph is 
endowed with a binary variable Sj = ±1. At each time step an agent i is selected at 
random along with one of its neighbours j and the selected agent takes the opinion 
of the neighbour, i.e Sj = Sj. In the case of a voter model on a d-dimensional 
hypercubic lattice, the spin-flip rate for the site x is given by[|] 


W^{s) 


1 



Sx(t) / ^ '5y(^) 
y^N{x) 


( 2 . 1 ) 


where s = denotes the state of the system at time t, Sx is the value of the 

spin on site x, N{x) is the set of its neighbouring sites, and r dehnes the timescale of 
the process. This particular form of spin-flip rate, which is just a constant times the 
fraction of disagreeing neighbouring sites, dehnes the so-called linear voter model. It 
is possible to dehne other voter-like models in which the spin-hip rate is not simply 
a linear function of the local ehective held hx = J2y&j\f{x) satisfy the 


Z 2 symmetry and have similar properties [2S11]- We will focus on the model with 
spin-hip rate (2.1) here. Note also that we are taking x G T/'. We will assume the 
lattice to be inhnite in all the calculations appearing further in this section, though 
we will be especially concerned with hnite size ehects in Sec. 

Equation (2.1) implies that this spin model has no bulk noise, he., if a site 
‘agrees’ with all its nearest-neighbours, its spin-hip rate vanishes. In this sense, the 
dynamics are similar to the zero-temperature Glauber ones. The consequence is 
that the ‘consensus’ states, he., the states in which all sites have the same opinion. 


f Frachebourg and Krapivsky [8] define the spin-flip rate with a factor 1/r in front of the 
parenthesis, they take r = d/d, and they therefore have an overall factor d/4. This coincides 
with our definition of Wx since we have a factor l/(2r), we choose rd = 2 and we also have a 
factor d/4 overall. There is, however, a difference with the choice made by Ben-Naim et al. who 
used a r that is half ours in their calculations [24] . 
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are ‘absorbing’ states. Indeed, if the system reaches one of the two consensus states, 
it will never leave it. 

However, this does not mean that the asymptotic steady state must be one of 
full consensus. In fact, it turns out that the coarsening process is not always effective 
in bringing the system towards a single-domain state, and whether it does or not 
depends on the dimensionality of the lattice. For d < 2 the system coarsens until 
ultimately reaching a single domain state, while for d > 2 there is an inhnite family 
of non-completely-ordered steady states mm- The discrepancy in the asymptotic 
regime reached above and below d = 2 will be further discussed in this section. 

The probability distribution satishes the master equation 

^P(s, t) = Y, t) - W^{s)P{s, t)) , (2.2) 


where s® is the conhguration that differs from s only in that the spin on the site x 
is reversed. One can then derive a set of differential equations for the n-spin time- 
dependent correlation functions {sx^-.-Sx^) = Xls'Sa:i---Sa:„P(s,t) andfindthat, since 
the update rule Wx is simply linear in the local spin, the equations for the correlation 
functions of different order decouple. 

The single-body correlation function or average magnetisation satishes [SI El] 

= -2(SxWx(s)} = ^Ax(Sx} , (2.3) 

clt 2ra 

where A^, denotes the discrete Laplace operator, 

d 

Ax{Sx) = —2d{Sx) + ^(Sjc+ei) + (Sx-eJ^ , (2-4) 

i=l 

and are the set of unit vectors dehning the lattice. In the inhnite system 

size limit or for periodic boundary conditions all sites satisfy this same equation. 
For hnite size systems with open boundary conditions the sites at the edges should 
be considered separately. 

Similarly, for the two-body correlation function one has [HI El] 


— 2^S3.'Sy ^H4;(s)-|-IFy(s)^ ^ ^ ^Ax N Ay^ (S x S y) .(2.5) 

Interestingly enough, Eq. (2.2) is mathematically equivalent to the master 
equation for a continuous-time symmetric random walk on Z‘^ with jumping 
rate As a result, the mean magnetisation per site, dehned as m{t) = 

lim 2 ,_,.oo L}'i{^x(t)) plays the role of the total probability for the walker 

and is thus a conserved quantity. The same result can be obtained by summing 
both sides of Eq. (2.3) over all lattice sites. Notice that while the magnetisation of 
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a specific system does change in a single update event, the average over all sites and 
over all trajectories of the dynamics is conserved. 

Consider a hnite system with an initial fraction p of voters in the +1 state and 
1 — p in the —1 state, so that the initial magnetisation density is mo = 2p — 1. 
Suppose that the system reaches consensus in which the state of magnetisation 
m = +1 occurs with probability p+i(p) and the state with m = —1 with probability 
1 — p+i(p). Then since mo = moo one has 2p — 1 = 2p+i(p) — 1, hence p+i(p) = p. 


Concerning again Eq. (2.3), by using the discrete Fourier transform of (sa,), one 


can prove that its general solution on an inhnite size lattice has the form |H], 


{Sx{t)) G T ^ ^ CTyNx-y \ 


( 2 . 6 ) 


where ay = {sy{t = 0)) and Jx is a shorthand notation for the multi-index modihed 


Bessel functions, Jx{u) = ^ with Xq, the usual modihed Bessel function 

of order a. 

If the initial conhguration is such that a single -|-1 voter sits at the origin and 
is surrounded by a “sea” of undecided voters (he. Sa;(0) = ±1 with probability 


1/2, for all X ^ 0) then, since ax = Sx,o, the solution to Eq. (2.3) reduces to 


(sxit)) = e By using now the asymptotic relation Xa{z) ~ C / ^/2nz, 

z 3> 1, valid for any real a, one hnds the asymptotic behaviour of the average site 
magnetisation, {sx{t)) ~ [27rt/(r(i)]“'^/^. Thus, a single voter relaxes to the average 
undecided opinion of the rest of the population. 

The last result is exact, but does not provide meaningful information on how 
the steady state of the system is reached. In this sense, a more interesting quantity 


is the two-body correlation function determined by Eq. (2.5). In order to solve this 


equation [8] one makes the assumption that at each time t the state of the system 
is translationally invariant, so that (sxSy) depends on the lattice vectors x and y 
only through their difference n = x — y. Then, by denoting Gn(t) = {sxit)sx+n(t)), 


Eq. (2.5) simplihes to 


— — AnGnit) 
at TCL 


(2.7) 


which should be solved subject to the boundary condition Go(t) = {s^(t)) = 1, for 
any t. In addition, it is natural to choose the initial condition G.ri(0) = that 
corresponds to a completely uncorrelated initial state. Equation (2.7) is basically 


identical to Eq. (2.3) apart from numerical factors, and one would be tempted to 


consider a solution of the form Gn{t) = e -r Nn {‘^t/{Td)). However, Go(t) does not 
satisfy the boundary condition. In order to maintain Go(t) = 1 throughout the 


evolution, one can reformulate the problem posed by Eq. (2.7) as the equivalent 
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lattice diffusion problem with a constant localised source at the origin, and look for 
a solution of the form 


Gn(t) — e 


2t I 2t\ 


dt' Sd{t — t') e 


"■(S) 


( 2 . 8 ) 


with Sd{t) the “strength” of the source. From a physical point of view, this solution 
corresponds to placing a source Go(t = 0) at the initial time at the origin and 
supplement it by an additional input Sd(t)dt that is added during the time interval 
{t,t + dt) to keep the overall value unchanged. Equation (2.8) evaluated at the 
origin (n = 0) becomes 


1 = 






^ 2,t 

Xn I —- 
Td 


(2.9) 


By using now the Laplace transform of the strength, Sd{X) = Jo dt Sd(t)e and 


the Laplace transform Td{\) of the function Td{t) = [Xo(t)e one arrives at 


S.(A) = -l + ^ 


x-taEx 


1 -1 


( 2 . 10 ) 


Using now the integral representation of the modihed Bessel function Xq, namely 
Mx) = ^ Jgit is possible to express Td in terms of the Watson integrals. 


r-27r 


r-27r 


u\) = 


(2vr) 


dgi ■ ■ ■ dqd 


X d 


( 2 . 11 ) 


22i=i cos qi 

and hnd an expression for Sd{X). For example, in the case d = 1, Si{X) = 
xJ{X + 2)/X. More complicated expressions arise when d is larger and ultimately 
there is no closed-form for them. Nevertheless, we are just interested in the 
asymptotic behaviour of the source strength Sd, which in turn is given by the low-A 
limit of its Laplace transform |H], 


Sd{X) 


and thus 


Sd{t) 


(-A) ^ 

(TA)"Mn“' |1 /(tA)] 

yN 


¥)■ 


if <i= 1 
if d = 2 
if d > 2 


as A —y 0 


( 2 . 12 ) 


In-(9 


if d = 1 

if d = 2 as t —)■ -l-cxD (2-13) 

const. if d > 2 

In d = 2, the long-time behaviour of the source strength in the integral is 
S' 2 (t —F) — l/ln[(t —F)/r] ~ l/ln(t/r). Using the asymptotic relations for X^,, and 
calling n = (^ 1 , 712 ), Eq. (2.8) implies 


i + taWT) 


dt' e' 


2t' 


t' 


T 


x„ J - 1 o I - 


■712 


t' 


(2.14) 
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as t —?• +C )0 dropping corrections with c a numerical factor to be determined. 

Using the integral representation of the modihed Bessel function X„ for integer values 
of n, that is X„(t) = [^cos/c — in k], Eq. (2.14) reduces to 

c 1 


Gr^(^) ~ 


dk^ 


dk 2 e 


—ink 


f{kO) + 0[j] (2.15) 


(2.16) 


ln(t/r) (27r)2 j 

where fc = (fci, k 2 ) and the function / (fc, t) is given by 

1 —-(2 —COSfci—008^2) 

/(fc,t) = r —- - -— . 

2 — cos ki — cos k 2 

Apart from a time-dependent prefactor, one can recognise in / the dynamical 
structure factor of the system, which is dehned as the lattice Fourier Transform 
of the space-time dependent correlation function, 

1 


S{kO) = 


f{kO) oc ^ G^(t) 


Ank 


(2.17) 


111 Wt) 

iLx:z^ 

In the limit |fc| —)■ 0, / can be approximated as f {k,t) — ‘2^ t k~‘^{l — e“^ ^^), where 
k = \k\, i.e. it becomes isotropic in fc-space. Then the large-distance behaviour of 
the correlation function is characterized by the scaling form 

^ ■(- 


GAt) 


f 


n 


1 ii («/ t ) \s/ij^^ 


(2.18) 


where the scaling function / is just given by the inverse Fourier Transform of /. 
Equation (2.18) clearly displays the emergence of a dynamical characteristic length 
i{t) which scales as -i/t, and the logarithmic violation of dynamic scaling. 

An interesting quantity that can be extracted from the two-body correlation 
function is the density of reactive interfaces p, dehned as the average value of the 
fraction of unsatished bonds or, equivalently, the fraction of neighbouring voters with 
disagreeing opinions. This quantity is linked to the correlation function through the 
relation 

- d 


Pit) = 2 




i=l 


= (2.19) 


where e, are the lattice unit vectors. Note that the sum over the nearest-neighbours 
can be lifted since the dynamics is isotropic along the d principal directions of the 
lattice. From Eq. (2.8) evaluated at n = (1,0, ...,0) and the fact that Go = 1, one 
obtains 


pd) = 2® 




\Td 


2t 


rd 


ZA 


2t 


1 

2 


du Sd{t — u) e ^ Xg ^ 


2u 

rd 


rd 


MS) 


MS 


( 2 . 20 ) 
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Combining the latter equation with Eqs. (2.13) and the asymptotic relations Xq 
Ii{z) ~ and Xo(z) — Xi( 2 ;) ~ I/x/Sttz^ -f 0{z~^/‘^) 

asymptotic behaviour of the density of reactive interfaces is found to be 


z) ~ 
, the 


\ ^ - 

p{t) ~ In ^(f/r) 
a — 


a d = 1 

ii d = 2 as t —)■ +oo ( 2 - 21 ) 

if d > 2 


These results allow us to establish some conclusions on the coarsening process 
in the voter model: in d < 2 the probability that two voters at a given separation 
had opposite opinion vanishes asymptotically, no matter how much distant they are, 
and coarsening eventually leads to a single-domain hnal state. In d > 2, an inhnite 
system reaches a dynamic frustrated state, where opposite-opinion voters coexist 
and continually evolve in such a way that the average concentration of each type of 
voters remains hxed. Dimension d = 2 is particular since it lies at the border between 
the two cases. There is a coarsening process which brings the system towards the 
single-domain state, but it is very slow, since the density of active interfaces vanishes 
only as l/ln(f/r). 

As a last effort, we derive the two-time correlation function, dehned as 
to) = {sx{t)sxito)), which is an interesting quantity to look at since it provides 
information on the typical timescale for the process to reach a steady state. For 
hxed to and Xq G let us introduce the function Fx(t;xo,to) = {sx(t + to)sxo(to)) 
for any a; G and t > 0. Dropping for a moment the dependence of F on 
to and Xq, it is easy to see that it satishes the same equation as the single-body 
correlation function, i.e. ^Fx{t) = ^^A^F^it), apart from a factor 1/2. Thus 
Fx(t) = e~^^'"Y.yfyJx-y{t/Td) where fy{xo,to) = {sy(to)sxo(to)). Then assuming 
that at each time the state of the system is spatially translational invariant and 
using Axit, to) = Exit- to; x, to), one gets 

Ax{t,to) = GM (^) > ( 2 - 22 ) 

n£Z‘i ^ 2 


with the dependence on x disappearing consistently with the hypothesis of 
translational invariance. As a simple check we verify that setting t = to in ( 2 . 22 ) 
we hnd Go(to)- Indeed, using J7n(0) = 11^=1 = 0 all n 7 ^ 0 and 

1 X 0 ( 0 ) = nf=i Xxi=o( 0 ) = 1 this fact is verihed. 

In the particular case to = 0, if the initial conhguration is completely 
uncorrelated, i.e. 0 ^( 0 ) = Snfi, the solution reduces to 

1 d 


Aoit) = A{t,to = 0 ) = e 




(2.23) 
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with asymptotic behaviour Ao(^) ~ [Int/. 

In the limit t 3> to one can use the asymptotic expansion of Nn{u) = 
Y[i=i^xX'^) — [e“/\/27rn]'^ with u = {t — to)/{Td) and, therefore, 

hmA3,(t,to) = [27r(t - to)/(Tci)]“'^''^ ^ Gn{to) . 


(2.24) 


The to-dependent last factor can be estimated as follows 


^{to) = Gn(to) •—t J 


d^x C{x,to) =drrC{r,to) (2.25) 


with C{x,t) the space-time correlation function in the continuum space limit, 
and using tl 

K{to) = 271 


Setting d = 2 and using the scaling function for G(r, to) expressed in Eq. (2.18) 

271 


dr 


ln(to/r) ^ ln(to/r) r ‘ 


to 


Going back to Eq. (2.24) this implies 

lim A^(t, to) oc , . 

ln(to/r) 


(t/to - 1) 


-1 


(2.26) 


(2.27) 


Further details on how to obtain the analytical results sketched in this section can 
be found in [HI ETj . 

We have already explained how the asymptotic behaviour of the space-time 
dependent correlation functions can be obtained in a way that exploits the special 
properties of the master equation. An alternative treatment of the many-body 
correlation functions uses an equivalence between the voter model and an auxiliary 
process of annihilating random walks [HHliniE]. By using this approach, Scheucher 
and Spohn obtained the same result for the dynamical structure factor in the small 
k and long-time limits for d=2 


S{k,t) 


ln(f/r) K^{t/T) \K{t/T) 


\k\ 


(2.28) 


with K{t) (X t 2 and F{u) = ^(1 — e “^), as found by employing the master 
equation formalism. From here one recovers the asymptotic form for C{x,t) in 
Eq. ([2T^. 


3. Numerical analysis 

In this Section we present our numerical results. We hrst compare them to the 
analytical ones recalled in Sec. for inhnite size systems and we later focus our 
attention on hnite size effects. 
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We define the model on a square lattice with linear size L and periodic boundary 
conditions. In all cases we start the dynamics at time t = 0 with a random initial 
condition with Sa, = ±1 with probability a half. 

One unit of time {i.e. r) corresponds to spin-flip attempts. As the system 
coarsens the number of flippable spins decreases and more attempts are necessary 
to change the conhguration signihcantly. In order to accelerate the simulations we 
used a continuous time Monte Carlo algorithm with the voter model dynamic rule. 
Unless otherwise stated, the quantities that we present below were averaged over 
10^ samples. 

As we will be particularly concerned with the geometric properties of the 
coarsening process, let us give here a number of dehnitions that we will use in 
the rest of this Section. We dehne a cluster or geometric domain as the ensemble of 
hrst-neighbour parallel spins. The cluster area is the number of spins belonging to 
it. Any such domain is surrounded by an interface that corresponds to the ensemble 
of broken hrst-neighbour links surrounding the cluster. The total interface length 
(external plus internal) is the number of such oppositely oriented spin pairs. 

3.1. Snapshots 

In Fig. [^we show three series of snapshots of the bidimensional voter model (hrst 
row) and the ferromagnetic Ising model at times f = 4, 64, 512, 4096 (henceforth 
all numerical times are expressed in MCs and we omit this time unit to lighten the 
notation). The Ising model (IM) has been quenched to zero temperature (second 
row) and the critical point (third row) and it evolves with a heat-bath Monte 
Carlo algorithm. Red and white points represent the two spin conhgurations. The 
snapshots illustrate the coarsening phenomenon induced by the different microscopic 
dynamics. In the case of the 2d\M. instantaneously quenched to T = 0 the dynamics 
are purely curvature-driven: for sufficiently long time, all the interfaces move with 
a local velocity that is proportional to the local curvature [la m. As a result 
the interfaces tend to disappear independently of one another, i.e. there are no 
coalescence processes. Instead, in the voter model the dynamics are driven by 
interfacial noise. In other words, if the initial conhguration consisted of a single hat 
interface between two domains of opposite opinion, opinions would slowly dihuse 
from one domain into the other and, after a sufficiently long time, the original sharp 
interface would become a dihuse interface. As one can see from the snapshots, 
phase-ordering still occurs but the resulting domains are very jagged and preserve 
their fractal geometry even at the late stages of evolution. Note, however, that 
the dynamics of the zero temperature Ising model and the voter model have one 
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important feature in common, namely, they are both characterised by the absence of 
bulk fluctuations. But they also show one important difference in the morphological 
properties of their interfaces. Indeed, the domain walls in the voter model are more 
similar to the ones in the critically quenched Ising model, shown in the third series 
of snapshots in the same hgure, than to the ones in the Ising model evolving at any 
subcritical temperature. The critical conhgurations are, though, plagued with bulk 
fluctuations, and these are absent in the voter model. 



(a) t = 4 (b) t = 64 (c) t = 512 (d) t = 4096 



(e) t = 4 (f) t = 64 (g) t = 512 (h) t = 4096 



(i) t = 4 (j) t = 64 (k) t = 512 (1) t = 4096 


Figure 1: Snapshots of the voter (first row) and Ising (second and third rows) models on a 2d 
square lattice with linear size L = 640 and periodic boundary conditions. The Ising model has 
been quenched to T = 0 (second row) and E (third row). The images were taken at the times 
indicated below the snapshots. 


In Fig. 1^ we display a series of snapshots of the voter model for even longer 
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(a) t = 4096 (b) t = 8192 (c) t = 16384 (d) t = 32768 



(e) t = 65536 (f) t = 131072 (g) t = 262144 (h) t = 524288 


Figure 2: Snapshots of the voter model with linear size L = 640 and periodic boundary conditions. 
The times shown are ti = 2*, with i = 12 — 19. With the same convention as in Fig. [l| the +1 
voters are coloured in red, while the —1 in white. We highlight here the percolating clusters with 
different colours. Percolating clusters of +1 (—1) opinion are green (blue). A single domain of 
“+1” voters was reached at t « 6.7 • 10®. 


times than the ones nsed in Fig. highlighting the percolating clusters of the two 
types. These conhgurations could be compared to the ones shown in [23] for the 
2(iIM quenched to T = 0. We note that the identity and form of the percolating 
clusters are not preserved in the hrst seven snapshots, until the system enters the 
late stage of evolution and hnds full consensus. 

In Sec. H we identify the largest cluster in the system and we study several of 
its properties letting us obtain in this way the exponent linking the system size to 
the time needed to reach percolation, and the fractal dimension of the percolating 
cluster and the one of its perimeter. 

In the case of a hnite lattice with periodic boundary conditions one can 
distinguish two types of domains: the ones that are homotopic to a point on the 
torus and the ones that wrap around the hole and cannot be completely shrunk 
without breaking into disconnected pieces. Even though the former can have a 
linear size comparable or even bigger than the one of the system, we identify the 
percolating clusters as the ones that wrap around the torus hole only. 
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Figure 3: Concentration of active interfaces p{t) = (1 — {sx{t)sx+ei(t)))/2 as a function of time. 
Data for L = 160 and 640 are presented in a log-linear form. The error-bars are smaller than the 
symbol sizes. The dotted (black) line represents the analytic asymptotic form p{t) = a/(21nt -|- 6) 
with a = 3.14 and b = 5.54. In the inset we show b = 7r/p(t) — 21nt as a function of t for different 
lattice sizes together with the analytic value 5.54 shown with a horizontal dashed line. See the 
main text for a discussion. 


3.2. Interface density 


In Sec. 1^ we provided an expression for the long-time behaviour of the fraction 
of active interfaces. In the disordered initial condition p(0) ~ 1/2. In d = 1 
the voter model is equivalent to the Glauber IM and pit) decays as while 

in d > 2 the density of interfaces converges to a constant, pit) a-ht-^!’^. In 
d = 3 the model has, then, blocked conhgurations asymptotically, as in the 2d and 
3d IM at T = 0 ^7\ |2B]. In the case of d = 2, pit) decays logarithmically and 
several authors [29], [30] tried to study this particular behaviour with Monte Carlo 


simulations. By starting from Eq. (2.20) it is possible to obtain a more rehned 
estimate of pit) [8], 


pit) = 


TV ^ / In f 

2Inf -|- In256 \ t 


(3.1) 


This result has to be contrasted to the algebraic decay, pit) ~ of curvature 

driven domain growth. For instance, in the 2dIM model this same quantity decays 
as pit) ~ with lit) the characteristic growing length and 

Zd = 2 . 

In Fig. [^we present numerical data for pit) in a voter model with linear size 
L = 160 and 640 with times reaching t = 10®. The analytic result in the asymptotic 
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L T 


(a) 


(b) 


Figure 4: (a) The average consensus time (T) (red diamonds) for different values of the lattice 
linear size L in a log-log plot. The value of (T) has been computed over approximately 10® 
samples for the smaller sizes, and over about 7 • 10® samples for L — 640. We report also the 
approximate width of the distribution, computed as the standard deviation of the data collected, 
ST=[N-^ as vertical dashes with length 2ST on each data point. Roughly, 

6T/(T) — 0.74 for all the sizes. The function C was fitted to the few data points available 
yielding 77 ~ 2.21 and C ~ 0.75. By introducing a logarithmic prefactor, Cx'^lnx, the fit yields 
instead 77 ~ 2.05 and C ~ 0.33. In the upper part of the same panel we report the skewness 7 and 
the kurtosis n of the sampled distributions in order to quantify the deviations from a Gaussian 
form, (b) The distribution of consensus times for different system sizes L given in the key. 


limit f{t) = a/(21nt -|- h) with a = vr and b = In 256 ~ 5.54 accompanies the data 
as a dotted (black) line. We have performed detailed fits of the data finding that 
the parameter a approaches the analytic value quickly. We then hxed a = n and 
we measured the parameter b by studying b = 7i/p{t) — 2 In t as a function of t for 
different system sizes {p{t) are the measured values). We show the result of this 
analysis in the inset to the same hgure. The approach to the analytic value shown 
as a black dotted horizontal line is indeed very slow. This fact explains why several 
authors did not see the asymptotic law in their numerical data and used instead a 
different logarithmic form, p{t) = Cln“'^f, with an effective exponent a ~ 0.6, to 
£t their data [SI EHl El] • 

3.3. Consensus time 

With much longer simulations we were able to measure the average consensus time, 
(T) = N~^ 'Ti with Ti the time required by the Tth sample to reach full 
alignment and Ng the total number of samples. In Fig. (a) we show the results 
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obtained by averaging over at least 7 ■ 10^ realisations of the dynamics for each value 
of the lattice size L. The averaged consensus time approximately follows the law 
(r)(T) However, in Sec. [^we recalled that the correlation functions suffer 

from logarithmic correction. Therefore, we tried to take into account this kind of 
correction by htting the function C In L to the data {L, {T){L)) and we obtained 
rj = 2.05 ± 0.01, C = 0.33 ± 0.02, and a better agreement with the numerical data 
than with the pure power law. 

An estimate of the characteristic width of the probability distribution of the 
consensus time is given by the standard deviation, 6T = O'i ~ 

The relative standard deviation 5T / (T) was found to he in the interval 0.70 — 0.75 
for every L and no particular dependence on the number of samples was observed 
for Ng > 10^. To highlight this behaviour we added vertical dashes of width 26T{L) 
centred on each one of the data points in Fig. (a). We stress that these dashes do 
not represent any type of error on the numerical value of the average consensus time, 
but only a measure of the average dispersion of our data on the available population 
of samples. 



r/{L" luL) 


Figure 5: Rescaled histogram of the consensus time pct{T, L) for different lattice sizes. Time T 
is rescaled by InL, while the value of pct{T,L) is multiplied by L“. The best collapse of the 
different curves was obtained for a ~ 2.22. 


In Fig. 1^ (b) we show the histogram of consensus times PcriT, L) for the 
different lattice sizes that have been simulated. The curves have approximately 
all the same shape when plotted in a log-log scale, so it is reasonable to assume the 
following scaling Ansatz 

T 


PcT{r,L) = L-^V 


LP InL 


(3.2) 
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with exponents a and /3 to be determined. Equation (3.2) simply states that the 
probability distributions of the consensus time for finite systems with different size 
have identical form, apart from a prefactor, if time is rescaled by a typical time 
Ttyp ~ InL. This is a natural choice since (T) ~ \nL and thus we expect 
(3 to be very close to 2. In Fig. [^we show the result of the scaling for (3 = 2 and 
a = 2.22. We could not appreciate signihcant improvements of the collapse for 
values of (3 slightly different from 2 , so our assumption was conhrmed. 

The upper part of Fig. (a) displays the skewness, 7 = Ng^J2f=i{'7i — 
(T))^/(5T^, and kurtosis, k = ~ {T)Y/5T^ — 3, of pcTiT^L) as a 

function of system size. The deviation from zero of the former quantihes the 
asymmetry of pcT and of the latter gives a further idea of the non-Gaussian 
character of the distribution. The skewness seems to have converged to a system 
size independent value that is slightly higher than one, while the kurtosis is still 
varying signihcantly for different system sizes. The last data points, for L = 640, 
are clearly not converged and many more samples would be needed to reach a good 
estimate for them. 


3.4- Persistence and autocorrelation 

In general one dehnes Pn{t) as the probability distribution for the number of 
opinion changes n experienced by a voter during the time interval ( 0 ,t), with 
n G {0,1,2,...} 121]. The hrst of these quantities, Poit), is equal to the fraction of 
voters who did not change opinion up to time t, i.e. the persistence probability [12|- 

In terms of spins, it measures the fraction of sites that have not experienced any 
spin-flip up to time t. In most statistical physics models ra. the persistence decays 
in time with a power law t~^ with a new independent persistence exponent 9. In the 
two-dimensional Glauber-Ising model at zero temperature, the exponent 6 has been 
evaluated numerically with high precision and it takes the value 6 ^ 0.199(2) for 
initial conditions with short-range correlations [33]. The asymptotic behaviour of 
Pq in the voter model in d > 2 was hrst found numerically 12a and then computed 
analytically with a mapping onto a continuum reaction-diffusion process and the 
use of held theoretical tools [31]. In d = 2, 

Po{t) ^ k ■ exp [—a\iY t + O {lnt)~\ (3.3) 

for our choice rd = 2. The diherence in the behaviour of the persistence between the 
2dIM and the 2d voter model was investigated by Drouhe and Godreche [25] who 
introduced a class of stochastic processes on a 2 d square lattice that interpolate 
between these two. They also conhrmed the unusual time-dependence of the 
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Figure 6: (a) The logarithm of the persistence probability, InPo; against \rT t, for a system of 
size L = 640. A linear ht to the data in this scale, f{x) = b — ax, gives a ~ 0.26 and b ~ 0.36 (the 
fitting curve is showed as a dashed blue line). The persistence has been calculated as an average 
over 5 • 10^ realisations of the dynamics. Data for L = 160 are shown in the inset along with 
their fitting curve, (b) Time-dependence of the autocorrelation with the initial configuration, Aq, 
in a double logarithmic scale. The function f{t) = at~^ has been htted to the numerical data, 
yielding b = 0.995 and a = 0.308, which are both in agreement with what theory predicts, i.e. 
A{t) ^ (7rt)~ , shown with a dark dashed segment. In the inset we show data for L = 160 analysed 
in a similar fashion. 


persistence decay in the 2d voter model with numerical simulations. 

In this context, we tried to recover the theoretical prediction in Eq. (3.3) with 
simulations of the voter model with different sizes. We present data for L = 160 and 
L = 640 in Fig. (a). By htting the simulation data to the function in Eq. ( |3.3[ ) we 
found a ~ 0.26 and k ^ ~ 1.44. The estimated value of a is quite close to the 

theoretical value predicted by Howard and Godreche lai, who found a ~ 1/4 with 
corrections of order 0.01. 


In Sec. we showed that the autocorrelation with a completely uncorrelated 
initial conhguration, Ao{t), has the asymptotic behaviour Ho(f) — in d = 2, 

with the choice r = 1. In Fig. we show numerical data for Ho(f) in a system 
with linear size L = 640. As one can see, the data are in good agreement with the 
theoretical prediction. 

In Fig. we plot, instead, the two-time autocorrelation function A{t,to) for 
values of to > 0, and two different lattice sizes. At sufficiently long waiting-time 
to the curves tend to flatten, losing their decay. This is clearer in the left panel 
where data for the smaller size, L = 160, are shown. In this case, all curves reach a 
plateau for t — to > 10^, signalling that the steady state has been reached. Indeed, 
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(a) 


(b) 


Figure 7: The two-time autocorrelation function A{t,to) against time-delay t — to for different 
values of the waiting-time to. Panels (a) and (b) show data for systems with linear size L = 160 
and L — 640, respectively. 


we have calculated the average consensus time for a system with size L = 160 (see 
Fig-i, and we found Tieo — 6 ■ 10^, which is compatible with the behaviour of the 
autocorrelation function. The same feature is expected to arise for the larger size 
L = 640 at a still longer time delay. In the case L = 640 we scaled the data to the 
analytic form (2.27) by plotting Into against t/to — 1 in Fig. The scaling 

is very good for to > 256. 



t/to - 1 


Figure 8: The two-time correlation function A{t,to) times Into against {t/to ~ l)j for different 
values of the waiting time to, in a log-log scale. A segment with slope —1 is shown as a guide-to- 
the-eye to stress the good agreement with the exponent found analytically. 
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The numerical analysis of several averaged correlation functions that we have 
presented so far is in good agreement with the theoretical predictions for inhnite 
size systems recalled in Sec. However, we will see in the following part of this 
Section that by studying other geometrical observables, we get access to aspects of 
the dynamics that remain hidden in the correlation functions. This analysis will 
allow us to uncover another dynamic regime. 

3.5. Averaged number of wrapping clusters 

We analysed the time dependence of the average number of wrapping domains per 
sample, iVp(f; L), by supposing that it only depends upon a scaling variable, 

NAt-,L)=N{J^y (3.4) 

with ap a. parameter to be determined and A/'(m) a scaling function. At strictly 
zero argument A/'(0) = 0 since the initial random conhguration is below the critical 
percolation threshold. At inhnite value of the argument A/'(oo) —)■ 1 since the hnal 
state of full consensus has a single domain. In order to estimate ap we tried to 
collapse Np{t]L) for different values of L by plotting it against the rescaled time 
t/L'^p with trial values of ap. The values of ap that gave the best collapse were 
found in the range 1.65 — 1.68, and in Fig. |^we present the case 

ap = 1.667 . (3.5) 

Even though the quality of the collapse for < 0.1 is very good, we must point 

out that this is a quite rough estimate, since we could not appreciate remarkable 
differences between slightly different values of ap in the aforementioned interval. 
Deviations from the desired scaling form Np{t]L) ~ are observed for 

^0.1. These are indeed expected since the system enters the next dynamic 
regime of approach to full consensus. As u increases from zero, f/{u) increases 
monotonically up to a certain value greater than 1. At this stage of the dynamics 
there are, then, states with more than one wrapping cluster. The scaling function 
next decreases converging to 1 from above. The exponent ap sets the typical time 
required for the system to reach a regime with wrapping clusters to tp ~ 


3.6. The largest cluster 

We identihed the largest cluster at each step of evolution and we computed several 
of its properties. This analysis allows us to distinguish whether the largest cluster 
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Figure 9: Average number of wrapping domains Np{t,L) against rescaled time t/L°'^, with 
ap = 1.667, for systems of sizes L = 20, 40, 80, 160. As t —>■ +oo the curves should converge to 
1 (single domain state). 


has wrapped around the system in one or more directions and, moreover, to which 
kind of criticality it belongs. 

We hrst measured the averaged number of its interfaces with positive, vanishing 
and negative curvature, (iV+), (Nq) and (iV_), respectively. A non-percolating 
cluster has a single external interface with positive curvature; we therefore call 
(Nep) = (Np) (with ep for external perimeter). A cluster that percolates in one 
direction has two interfaces with vanishing curvature. The interfaces with negative 
curvature are internal to the cluster and surround its holes. 

Another interesting observable is the area of the largest cluster Ac that we 
normalise as Ac/L^^ with Da a fractal dimension that we need to hnd. We recall 
that the fractal dimension of cluster areas in 2d site percolation [35] is 

~ 1.896 . (3.6) 

Finally, we calculated the averaged total length of the boundary C as the sum 
of the length of external and internal interfaces described above. We also normalised 
this length as with Dh a fractal dimension. For the sake of comparison, we 

recall that for 2d site percolation the cluster hull fractal dimension [36l ET] takes the 
values 

= 1.750 . (3.7) 

All these quantities show different scaling properties at short times before 
percolation is reached and at longer times, when the percolating cluster has 
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established, and the further evolution takes the system to the dual absorbing state. 
In the former time-regime, the data scale as a function of with ap = 1.667, 

see all panels in the right column in Fig. while in the latter they do as a function 
of (apart from logarithmic corrections), see all panels in the left column in the 
same hgure. The value of ap is rather large in the voter model (much larger than in 
the 2dIM where ap = 1/2 [23]) and this makes the distinction between the regime of 
approach to percolation and the further evolution towards full consensus very hard. 

The averaged number of interfaces with positive curvature smoothly decays 
from one to zero as more and more samples wrap around the sample in at least 
one direction, see Fig. 10 (a)-(b). Concomitantly, the averaged number of maximal 
size clusters that wrap around the sample in one direction increases in time from 
zero initially to a value that is close to 0.9 to later decay again to zero when the 
cluster percolates in the other direction as well, leaving only internal interfaces in the 
system (not shown). The averaged number of internal boundaries has a very similar 
qualitative behaviour to the one of the zero-curvature ones (not shown either). Panel 
(a) confirms the scaling with at long times while panel (b) indicates that the 
good scaling variable at short times is 

Figures [l0| (c)-(d) and (e)-(f) display the area of the maximum size cluster and 
the total perimeter length, respectively. The asymptotic approach of Ac /in (c) 
to one and C in (e) to zero conhrm that the systems approach full consensus. They 
also prove that there are no blocked states in the voter model as there are in the 


2dIM atT = 0[3Hl|39lIlQlllIlll2]. 

The scaling of the horizontal axis in panels (c) and (e) is intended to work at 
long times only. The persistent deviations from scaling should be due to the fact that 
we did not include in the scaling variable corrections that depend logarithmically 
on the system size, although we know from the analysis of the consensus time that 
these should exist. 

The scaling of the vertical axis in panel (d) is intended to determine the fractal 
dimension of the area of the largest cluster, in the short time regime. We used 
different values of Da and we found that the most satisfactory collapse in the 
short-time regime is found for Da in the interval 1.89 — 1.93. In panel (d) we 
show the scaling found using the fractal dimension of 2d critical site percolation. 
Da = D^^ ~ 1.89583 which is within this interval (while Da for other critical states 
lies outside this interval, i.e. Da = 1.948 at the critical Ising point). Once again, due 
to the large value of ap the two relevant asymptotic regimes are not well separated, 
and we cannot do better than this in the determination of Da- 

Similarly, the scaling of the vertical axis in panel (f) should determine the fractal 
dimension of the boundary of the largest cluster. Dp (although we stress that we 
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show data for the total interface length here). Using Dh = = 1.75 we see 

that all curves cross at t/L^^ ~ 1 suggesting that critical percolation is reached at 
around fp ~ Other choices for the value of Dh do not allow for data collapse 
at any value of the scaling variable. Furthermore, one could argue that the scaled 
data for t/L°‘^ are slowly approaching, for increasing system size, a flat form. 

In the inset to panel (f) we show, for comparison, the same scaling plot, 0/ 
against for the zero temperature Ising model quenched from Tq —)■ cxd, with 

Dh = 1.75 and ap = 0.5. Apart from the rather small system sizes, L = 40, 80, the 
data for the larger systems show a very good collapse over a rather large time- 
window. In the Ising model the two asymptotic time regimes are rather well- 
separated (as ap = 0.5 is quite different from za = 2) and this fact contributes 
to the good data collapse. 

3.7. Two growing lengths 

We conclude that, as in the 2d Ising model with non-conserved order parameter [2111 
EHESlEs], the dynamics of the hnite size 2d voter model takes place in two distinct 
regimes: the system hrst develops wrapping clusters of critical site percolation kind; 
once these are established, the further growth is a more usual coarsening process. 
The two growing lengths controlling the evolution in the two regimes are 

ip{t) ~ , icit) ~ . (3.8) 

In the Ising model on a square lattice the exponents ap and Zd are rather different, 
ap = 1/2 [23] and Zd = 2. Since ap is so small, for the system sizes used in numerical 
simulations the approach to percolation occurs in a few steps and the time window 
is not sufficiently long to allow for a careful dynamic scaling analysis. Instead, in 
the voter model, ap ~ 1.667 is quite large and not very different from Zd = 2. The 
system takes much longer to reach critical percolation and the advantage is that a 
rather wide time window can be explored in which the relevant growing length is 
ip{t). 


3.8. Density of domain areas 

We now turn to the statistics of domain areas. We recall that we dehned a cluster 
or a domain as the connected ensemble of nearest-neighbour parallel spins and its 
area as the number of spins in it. 

Figures [TT] shows raw data for the time-dependent number density of domain 
areas in systems with linear size L = 640. Initially, the curves show no particular 
structure, as the random initial condition is non-critical and the weight of the 
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distribution at large areas drops significantly. However, as time elapses, a power 
law extending over several decades develops, as already noticed by Scheucher and 
Spohn m- Moreover, a bump with support over areas that are of the order of 
magnitude of the size of the system also appears. These areas are the ones of clusters 
that wrap around the system, as the ones discussed in the previous subsection. The 
height of the bump increases in time. It tends to become stable at the longest time- 
scales used, t ~ 17 ■ 10^. (For a smaller system size, say L = 160, the same features 
are realised but, in contrast, after growing in height the bump tends to wash out, 
after times of the order of 2 ■ 10^.) This feature is very similar to what was observed 
in the 2dlM quenched to zero temperature, although a stable bump, linked to the 
system reaching critical percolation, establishes at a much shorter time-scale, tp ~ 5 
for similar system sizes [201 EH as ap = 1/2 in this case [23] • Indeed, the curves for 
f = 0, 1, 4 in the inset are qualitatively identical to the ones for f = 0, 8,16384 
in the main plot. (We will make a quantitative comparison between the behaviour 
in the two models below.) In the case of the 2d voter model we observed that 
percolating clusters can appear in early stages of the dynamics, but they tend to 
break soon after their formation and reappear later on, taking longer to establish 
a stable pattern, see Fig. In particular, for a system with linear size L = 640, a 
stable structure of percolating domains establishes only after a time of the order of 
tp ~ lOF 

The analytic and numeric analysis of the 2dlM quenched from infinite to 
zero [201 EH or the critical m temperature showed that the number density of 
areas approaches a scaling form 


A'' nd{A,t) = d) 






(3.9) 


where Da is the fractal dimension of the areas studied, i{t) is the relevant growing 
length and <I> a proper scaling function. After tp this scaling has to be corrected by 
an additive term that takes into account the percolating clusters that had already 
established (the bump). In the zero temperature quenches of the 2dlM the approach 
to percolation was so fast that the study of this scaling for times such that the 
relevant growing length is ip{t) was not performed. In the critical quenches in [22] a 
triangular lattice for which the system was already at the critical percolation point 
initially was used. Here, with the voter model, we have the possibility of studying 
the dynamic scaling in the regime of slow approach to percolation in detail, by 
taking advantage of the large value of ap. 

The same datasets used in Fig. 11 are plotted in the form A^nd{A,t) against 
A in Fig. [^with u ~ 1.98 for L = 640. The value of the exponent u was found by 
fitting the data nd{A, t) corresponding to the longest time reached in the simulation 
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(t = 2048 for L = 160, not shown, and t = 16384 for L = 640) with the power 
law CdA~^, in the range of areas [10^, 10^] for L = 160 and [10^, 10^] for L = 640. 
We found Cd = 0.0207 ± 0.0001 and v = 1.972 zb 0.005 in the case L = 160, and 
Cd = 0.0220 ± 0.0001 and z/ = 1.980 ± 0.001 for L = 640. 

The value of the exponent v increases very weakly with L and should be larger 
than 2 in the inhnite size limit to ensure that the average area of non-percolating 
domains, / dA A nd{A,t), is a hnite quantity. However, the approach to the 
asymptotic limit is so slow that it is very hard to get closer to it numerically. This 
particular feature of the dynamics was also observed in HH for the voter model and 
we stress that, in the 2dIM, the expected value u = 2.05 is found only for a very 
careful choice of the areas to ht. 

The value taken by the constant Cd is very relevant to our discussion. Indeed, 
it was used in [20l [21] to distinguish the criticality of large scale domains in zero- 
temperature quenches of the 2dIM from equilibrium at Tq —?• oo and Tq = Tc. 
More precisely, the area distribution of clusters of occupied sites at critical site 
percolation and, say, domains of positive spins at the critical Ising point are given 
by nd{A) ~ 2CjA^^^ and nd(H, 0) ~ CIA^''\ respectively, with 2C ~ 0.023, and 
= 1 - 1 - d/D'^J and = 1 -b d/D^^ the Fisher exponents related to the fractal 
dimensions of the domain areas under the two critical conditions. Cardy and Ziff |13] 
obtained these universal constants analytically for the number density of hull- 
enclosed areas instead of domain areas using a Coulomb gas approach. Arguments 
presented in IZH suggest that very close values should apply to domain areas as 
well. Numerical simulations on the square and triangular lattice conhrmed the value 
obtained with held theoretic methods for hull-enclosed and domain areas HaiEilEs] 
After a zero-temperature quench of the 2(iIM with initial states drawn from inhnite 
temperature and critical temperature conditions, the evolving large scale areas 
are distributed algebraically and the number densities have Fisher exponents and 
constants in the numerator that are the ones cited above for critical site percolation 
(see the inset to Fig. 13) and critical Ising conditions, though both multiplied by a 


factor of two when clusters of both (up and down) species are counted 

In the voter model with random initial conditions we hnd Cd — 0.022 that is 
consistent with Cd = 2C (within numerical accuracy), see Fig. and with critical 
Ising initial conhgurations we hnd a constant taking the value q/2 = C (not shown). 
This result conhrms the reduction of the number density of hnite areas by two for 
initial conditions with long-range correlations with respect to the ones with only 
short-range correlations. 

In Fig. 13 we show A'' nd{A,t) against the rescaled area A/C for systems with 
L = 640. The value of a that allowed us to obtain the best collapse was found to be 
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approximately equal to 1.19. This value is in good agreement with the prediction 
ip^{t) ~ t^A/ap — ^1.90/1.67 _ ^1.14^ jnset we perform the same analysis 

on the 2dlM quenched to T = 0, by focusing on the very short time dynamics 
such that t tp. In agreement with the proposal, the curves collapse if one uses 

= ^/^( 187 / 96 )/ 0.5 ^ 

Apart from deviations caused by the appearance of wrapping domains in the 
bump, for large values of t all curves seem to have the same behaviour, namely 
two distinct regimes: for A/C ^0.1 there is a nearly flat region, which would 
mean that <!>(«) ~ const, and thus nd{A,t) ~ A~'^, i.e. the statistics of domain 
areas is independent of time; for A/C >0.1 instead, the scaling function seems 
to behave as an increasing power law <!>(«) ~ cm“, with a self-similar statistics of 
domain areas in the sense that it depends only on A/C. A £t of the data for the 


shortest time in Fig. 13 on the interval [1,100] of the scaling variable u = A/C, 
yields c = 0.034 ± 0.001 and a = 0.257 ± 0.001. Analogously, for the case L = 160 
we obtained c = 0.032 ± 0.001 and a = 0.264 ± 0.001 (not shown). The existence of 
two distinct regimes for small and large values of A/C is also observed in the Ising 
model as shown in the inset of Fig. Moreover, for large values of A/C we also 
observe an increasing power law (also shown in the inset) with a very similar power, 
i.e. c = 0.024 ± 0.001 and a = 0.262 ± 0.001 for L = 640. 

Coming back to the domain area statistics, as one can see from the plots, the 
flat region for A/C < 0.1 becomes larger as time increases, while the complementary 
region of larger domains shrinks, until disappearing. Note that as time increases the 
percolating (wrapping) domains become more and more predominant and eventually 
the number domain density converges to the absorbing state form which is just a 
delta function centered at A = Even though this fact does not rule out the 
possibility of a transient regime in which more than one stable percolating clusters 
coexist, similarly to what happens in the zero-temperature 2dIM. on a hnite lattice, 
we found that it establishes during a very short time period (compared to the whole 
duration of the dynamics) before the consensus state is reached, so it is somehow 
difficult to “catch” it in the domain area statistics. 


3.9. Space-time correlation function 

Having established the existence of two dynamic growing lengths in a hnite size 


system, we now put the scaling form of the space-time correlation, Eq. (2.18), to 


the numerical test. Figure 14 (a) shows data for C{r,t) on a lattice of linear size 
L = 640. The correlation function was calculated only along a principal direction of 
the lattice (e.g. the horizontal direction), as C{x,t) should be isotropic and depend 
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on X only at distances much longer than the lattice spacing. In Fig. 14 (b) the 
correlation function multiplied by Int is plotted against the scaled distance x/\/i. 
The curves at different times tend to collapse even though they deviate for large 
values of x/\/t. These deviations are due to hnite size effects: since we have taken 
periodic conditions at the boundaries, the data at distances of the order of the lattice 
size, specihcally x ~ L/2, are much affected by the boundaries. One reckons that, 
consistently, the deviation from the scaling law for large x/\/t occurs at smaller 
values of the scaling variable at longer times. 

However, from the analysis of the clusters we now known that in the dynamic 
regime in which percolating clusters develop there is another characteristic length 
in the problem, icit)- Accordingly, the scaling form of the correlation function has 
to be modihed to capture the dynamics in both dynamic regimes (see [23] for this 
analysis in the 2dlM). We therefore introduce a new two-variable scaling function 
g{u,v) 


C{x, t; L) = 


In it/r) 


\x\ 


ip{t/T) 


(3.10) 


Joit/r) ’ L 

such that for t tp L°‘^, the new scaling variable is close to one, ip{t)/L ~ 1, 
and one recovers the inhnite size limit. This suggests that the data for C{x,t] L) at 
different times ti and sizes Li chosen in a such a way that the ratio y = ip{ti)/Li 
is kept constant should collapse when plotted against the scaled distance x/lcit) ~ 
x/y/t, since Zd = 2 m. the voter model. In order to put this proposal to the test 
we computed the correlation function on square lattices with sizes Lj = 2*Lo for 
Lq = 80 and z = 0 — 4 and times ti = with to = 256, such that x = /Lq. 

As far as the exponent ap is concerned, we estimated it from the analysis of the 
largest cluster obtaining ap ^ 1.667, see Sec. 3.6[ and we conhrmed its value with 


the study of the time evolution of the number of percolating domains, see Sec. 3.5 


Figure 15 (a) and (b) presents the scaling forms in Eqs. (2.18) and (3.10), 
respectively. It is clear that the introduction of an extra scaling variable with the 
dependence on the new length scale allows us to achieve a much better data collapse. 


4. Conclusions 

The main goal of this work was to improve the understanding of coarsening in 
models with microscopic dynamics that are not driven by the minimisation of 
a thermodynamic potential and do not satisfy detailed balance. More precisely, 
we focused on a Z 2 symmetric lattice model with pairwise interactions driven by 
interfacial noise, viz. the 2d linear voter model on a square lattice. 
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We showed that the dynamic evolution of the bidimensional voter model on a 
square lattice proceeds in two distinct dynamic regimes. In the hrst one, the model 
approaches critical site percolation. The time needed to reach one such typical state 
diverges with the size of the system algebraically, tp — with the exponent 

ttp ~ 1.667 that is much larger than the one previously evaluated in the 2dlM 
quenched to T = 0 [22]. Next, the model evolves following a different mechanism in 
which consensus is progressively attained. The characteristic growing length for this 
process is also algebraic, i{t) ~ though with a different dynamic exponent, 

Zd = 2. In the social dynamics context, the hrst process can have important 
consequences. 

We based the conclusions above on the careful use of numerical methods. 
We hrst tested this approach against the theoretical predictions that were already 
available for inhnite size voter model. Most of the computed quantities, such as 
the fraction of active interfaces, the autocorrelation function and the persistence, 
were found to be in very good agreement with the analytic predictions for inhnite 
size systems. In particular, the peculiar logarithmic decay of the magnitude of the 
two-body correlation function and of the fraction of active interfaces was recovered, 
even though these results could be improved by simulating larger systems on very 
long times. We then focused on the spin conhgurations and from the analysis of 
their statistical and geometric properties we uncovered the approach to critical site 
percolation. Once the new growing length scale identihed, we used it to improve the 
scaling of the space-time correlation function for hnite size systems. 

In a series of papers, the role played by the approach to critical percolation in 
spin models with Ising [20l EU |22l ES] 01 EH] or Potts HHIET] variables in two or 
three dimensions jH] was studied. In all these cases the dynamics satisfy detailed 
balance and eventually take a hnite size system to thermal equilibrium. In this 
paper we explored a diherent kind of microscopic dynamics that does not satisfy 
detailed balance and approaches an absorbing state asymptotically. We still found 
a similar approach to critical percolation as in the ‘equilibrium’ cases albeit with a 
much slower growing length. 

In [23] we observed that, for the Ising model with microscopic dynamics 
satisfying detailed balance the exponent ap coincided with the ratio between the 
dynamic exponent for the late stage growth, Zd, and the lattice regular or averaged 
coordination number, Uc, i.e. ap = Zd/uc. In the voter model the dynamic exponent 
is Zd = 2 (on top, dynamic scaling for inhnite size systems suhers from logarithmic 
corrections) and, though a coordination number cannot be really identihed, we can 
claim that an ehective one is somehow larger than one. With this identihcation, the 
value of ap found numerically has the good trend, in the sense that the coordination 
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number is smaller than in Ising and this leads to a larger ap. 

This works opens several lines for future research. On the one hand, it would 
be interesting to extend this analysis to different types of lattices, variants of the 
update rule (with, e.g., local conservation laws [IH], memory [SU] or inhomogeneities 
in the form of zealots lai) and upgrading the voters to have many opinions (see 
and references therein). On the other hand, it should be possible to extract the 
growing length ip{t) analytically by taking into account the hnite size effects in the 
approach explained in Sec. 
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(e) (f) 


Figure 10: The largest cluster. The horizontal axes are t/Li^ and with ap = 1.667 in 

all panels in the left and right columns, respectively, (a)-(b) The averaged number of external 
interfaces (only for non-percolating clusters). Its area normalised by in (c) and with 
Da = 1.89583 in (d). (e)-(f) Its total perimeter length in the voter and 2dIM quenched to T = 0 
(inset) normalised by with Dp = 1.75. 























CONTENTS 


32 



Figure 11: Instantaneous domain area number densities. Main plot: 2d voter model with linear 
size L = 640. Inset: 2(iIM quenched toT — 0 from Tq —>■ oo and L = 640. In both cases the curves 
are presented in a double logarithmic scale and the times at which the data are collected are given 
in the key. 
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Figure 12: Time-dependent number density of domain areas, Ud, multiplied by A'' with z/ ~ 1.98 
for L = 640. As the system comes closer to the percolating state, the number density of domain 
areas tends to the time-independent form nd{A,t) Cd A'', save for a correction due to wrapping 
domains in the bump. The value of v was estimated by fitting the aforementioned functional form 
to the data corresponding to the latest time. The horizontal line is at Cd ~ 0.022. 
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^ 10-1 


10-2 


Figure 13: Number density of domain areas multiplied by A'' against the rescaled area Ajf^, for 
L = 640. The exponent v takes the same values as in Fig. v ~ 1.98. The value of a that yields 
the best time-scaling is a ~ 1.19. A fit to the function /(u) = c of the data at t = 64 over the 
interval [1,100] yields a « 0.26. The red dashed line was added to better visualise this power-law 
behaviour. Inset: the same scaling plot for the 2dlM after a quench to T = 0, at five short times 
given in the key. The plateau is at 2cd — 0.044 as explained in [5D]. The power law shown with a 
dashed line was drawn with the same value of the exponent a as for the voter model. 



(a) 
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Figure 14: In panel (a) log-linear plot of the space-time correlation function along a principal 
direction of the lattice, i.e. C{x,t) = {sxei{t)so{t)), at different times for a system with linear size 
L — 640. In panel (b) the same data as in panel (a) multiplied by Int and plotted against the 
scaled distance x/Ct in linear-log scale. See the main text for a discussion. 




Figure 15: C{x,t) ■ Int against scaled distance with Zd = “2 the dynamical exponent, 

in a linear-log plot. In panel (a) data are for L = 80 at different times given in the key. Panel 
(b) reports data for different Li and the corresponding times U such that jti is held constant. 
ap ~ 1.667 is the exponent in the relation tp ~ with tp the time needed for stable percolating 
domains to establish (see Sec. 3.6 for more details on tp). 



















